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ABSTRACT 


Time-varying dispersion and multipath propagation in a shallow underwater 
environment causes intersymbol interference in underwater communication. This thesis 
investigates a mitigation procedure for communication using a Binary Phase-Shift 
Keying (BPSK) signal. The method employed uses the time-reversed ocean impulse 
response to mitigate the degradation of the bit error rate performance. All results were 
achieved by the use of computer simulation of typical shallow water environments. 


v 



THIS PAGE INTENTIONALLY LEFT BLANK 


vi 



TABLE OF CONTENTS 


I. INTRODUCTION.1 

A. GENERAL.1 

B. OBJECTIVE.1 

C. ORGANIZATION.1 

n. BINARY PHASE-SHIFT KEYING.3 

A. BINARY PHASE-SHIFT KEYING.3 

B. SIGNAL CONSTELLATION FOR BPSK.5 

C. PROBABILITY OF BIT ERROR FOR BPSK.5 

D. BIT ERROR PROBABILITY CONVERGENCE.10 

E. TWO TYPES OF BIT ERROR PERFORMANCE DEGRADATION ....12 

HI. EXPERIMENTAL DESCRIPTION.15 

A. BPSK SIGNAL GENERATION.15 

B. MMPE MODEL DESCRIPTION.18 

1. A Brief Description of the MMPE Model.18 

2. Signal Representation in MMPE.25 

a. Time Windowing .27 

b. Processing of Bandpass Signals . 28 

3. Input Parameters for the MMPE Model.31 

C. OCEAN ENVIRONMENTAL CHARACTERIZATION.33 

1. Case 1: Positive SSP Gradient.33 

2. Case 2: Strong Negative SSP Gradient.34 

3. Case 3: Negative SSP Gradient Below Surface Duct.36 

D. BPSK DEMODULATION AND DETECTION.36 

1. Correlation Receiver..36 

E. MULTIPATH MITIGATION.39 

IV. SIMULATION RESULTS.43 

A. EVALUATION OF BIT ERROR PROBABILITY FOR BPSK 

SIGNAL.43 

1. Evaluation for BPSK Parameters.43 

2. Influence of AWGN on Bit Error Probability.48 

B. BIT ERROR DEGADATION AND MULTIPATH MITIGATION 

FOR BPSK SIGNAL IN A SHALLOW WATER ENVIRONMENT.50 

1. Bit Error Performance Results for a Positive SSP Gradient.50 

a. Investigation Of The One Specific Case . 50 

b. Results For Many Chosen Ranges and Depths . 58 

2. Bit Error Performance Results for a Strong Negative SSP 

Gradient.60 

3. Bit Error Performance Results for a Negative SSP Gradient 

below Surface Duct.62 

vii 





































V. CONCLUSIONS AND RECOMMENDATIONS.65 

A. CONCLUSIONS.65 

B. RECOMMENDATIONS.65 

APPENDIX A. MMPE INPUT FILES FOR THREE DIFFERENT OCEAN 

ENVIRONMENTAL CASES.67 

A. MMPE INPUT FILES FOR POSITIVE SSP GRADIENT.67 

1. pef iles . inp File of the Main Input File.67 

2. pesrc. inp File of the Source Data.67 

3 . pessp. inp File of the Environmental Data.68 

4. pebath. inp File of the Environmental Data.68 

5. pebotprop. inp File of the Environmental Data.69 

6. pedbath. inp File of the Environmental Data.69 

7. pef iles . inp File of the Environmental Data.69 

B. MMPE INPUT FILES FOR STRONG NEGATIVE SSP GRADIENT...70 

1. pef iles . inp File of the Main Input File.70 

2. pessp. inp File of the Environmental Data.70 

C. MMPE INPUT FILES FOR NEGATIVE SSP GRADIENT BELOW 

SURFACE DUCT.71 

1. pef iles . inp File of the Main Input File.71 

2. pessp. inp File of the Environmental Data.71 

APPENDIX B. DETAILED SIGNAL PROCESSING STEPS.73 

A. RANDOM BINARY DATA GENERATION.73 

B. FILTERING.74 

C. GENERATION OF BPSK SIGNAL.74 

D. DEMODULATION OF BPSK SIGNAL, FILTERING.75 

E. BER COUNTING.76 

APPENDIX C. COMPLETE PROCEDURES TO GENERATE THE PASSBAND 
OCEAN IMPULSE RESPONSE AND FREQUENCY RESPONSE FROM 
THE MMPE MODEL.79 

A. MODIFIED OCEAN FREQUENCY RESPONSE AND OCEAN 

IMPULSE RESPONSE FROM MMPE MODEL.79 

B. BASEBAND OCEAN IMPULSE RESPONSE BY PADDING ZEROS ..79 

C. FILTERING AND PRE-ENVELOPE OF THE OCEAN IMPULSE 

RESPONSE.81 

D. PASSBAND OCEAN IMPULSE RES PONSE.82 

LIST OF REFERENCES.83 

INITIAL DISTRIBUTION LIST.85 


viii 
































LIST OF FIGURES 


Figure 2.1. BPSK Signal in the Time Domain.3 

Figure 2.2. BPSK Signal in Frequency Domain.4 

Figure 2.3. Signal Constellation for BPSK.5 

Figure 2.4. Conditional Probability Density Functions: p(z\s 2 ) .8 

Figure 2.5. General Shape of the P B versus / N 0 Curve.11 

Figure 2.6. (a) Loss in E b /N 0 . (b) Irreducible P B Caused by Distortion.13 

Figure 3.1. BPSK Signal Generation.15 

Figure 3.2. Three Stages to Extract Ocean Response from MMPE Model.18 

Figure 3.3. Pulse Arriving after Propagation Delay x .26 

Figure 3.4. Spectrum of Bandpass Signal.27 

Figure 3.5 Spectrum of a Complex Valued Lowpass Signal Corresponding to 

Bandpass Signal of Figure 3.4.29 

Figure 3.6. pef iles . inp File of the Main Input File.31 

Figure 3.7. pesrc . inp File of the Source Data.32 

Figure 3.8. pessp . inp File of the Environmental Data.33 

Figure 3.9. (a) Positive SSP (b) Sound Propagation at Frequency 400 Hz, Source 5 m. ..34 

Figure 3.10. (a) Negative SSP (b) Sound Propagation at Frequency 400 Hz, Source 30 

m.35 

Figure 3.11. (a) SSP (b) Sound Propagation at Frequency 400 Hz, Source 5 m.36 

Figure 3.12. Two Basic Steps in the Demodulation/Detection of the Received Signal.37 

Figure 3.13. Correlator Receiver.38 

Figure 3.14 Multipath Mitigation for the Distorted Signal in Underwater Environment. ..40 

Figure 4.1. The Power Spectral Density for BPSK Signal.44 

Figure 4.2. General FIR Filter (Hamming Window Design).46 

Figure 4.3. Frequency Response for BPSK Signal Sampled by Null-to-null 

Bandwidth.47 

Figure 4.4. Frequency Response for BPSK Signal Sampled by Power Bandwidth.47 

Figure 4.5. BER versus SNR(dB).49 

Figure 4.6. BPSK Signal in Time Domain and Frequency Domain.51 

Figure 4.7. Passband Ocean Frequency Response (Magnitude and Phase) and Ocean 

Impulse Response.52 

Figure 4.8. Received Signal in Time Domain and Frequency Domain.53 

Figure 4.9. Comparison of Recovered Binary Data with Transmitted Binary Data for 

the Received Signal.54 

Figure 4.10. Time Reversed Ocean Frequency Response (Magnitude and Phase) and 

Ocean Impulse Response.55 

Figure 4.11. Mitigated Ocean Impulse Response and Frequency Response.56 

Figure 4.12. Mitigated Signal in Time and Frequency Domain.57 

Figure 4.13. Comparison of Recovered Binary Data with Transmitted Binary Data for 

Mitigated Signal Showing Zero Error.57 


































Figure 4.14. Bit Error Performance for Received Signal at Chosen Depths 5 to 95 m 

and Ranges 0.5 to 9.5 km for a Positive SSP Gradient.59 

Figure 4.15. Bit Error Performance for Mitigated Signal at Chosen Depths 5 to 95 m 

and Ranges 0.5 to 9.5 km for a Positive SSP Gradient.60 

Figure 4.16. Bit Error Performance for Received Signal at Chosen Depths 5 to 95 m 

and Ranges 0.5 to 9.5 km for a Strong Negative SSP Gradient.61 

Figure 4.17. Bit Error Performance for Mitigated Signal at Chosen Depths 5 to 95 m 

and Ranges 0.5 to 9.5 km for a Strong Negative SSP Gradient.62 

Figure 4.18. Bit Error Performance for Received Signal at Chosen Depths 5 to 95 m 
and Ranges 0.5 to 9.5 km for a Negative SSP Gradient below Surface 

Duct.63 

Figure 4.19. Bit Error Performance for Mitigated Signal at Chosen Depths 5 to 95 m 
and Ranges 0.5 to 9.5 km for a Negative SSP Gradient below Surface 

Duct.64 

Figure B-1. Random Binary Data in Time Domain and Frequency Domain.73 

Figure B-2. Frequency Response of the FIR Filter and the Filtered Binary Data.74 

Figure B-3. BPSK Signal s[n] in Time Domain and Frequency Domain.75 

Figure B-4. Demodulated BPSK Signal y[n] in Frequency Domain.76 

Figure B-5. Demodulated and Filtered BPSK Signal in Time and Frequency Domain.76 

Figure B-6. Comparison between Recovered Binary Data and Transmitted Binary 

Data.77 

Figure C-l. Modified Ocean Frequency Response H (/) and Impulse Response 

h[n\ .80 

Figure C-2. Baseband Ocean Frequency Response H 0 (f ) and Impulse Response 

A l«l.80 

Figure C-3. Frequency Response of FIR Filter and Filtered Baseband Ocean 

Frequency Response.81 

Figure C-4. Pre-Envelope Ocean Frequency Response H (f ).82 

Figure C-5. Passband Ocean Frequency Response H b (/).82 


x 



















LIST OF VARIABLES 


AWGN additive white Gaussian noise 

A c sinusoidal carrier amplitude 

{a n } message data or a set of random variables with a n = ±1 

BPSK binary phase-shift keying 

b(t ) random binary signal 

B T (/) Fourier transform of the truncated signal, b T (t) 

B(f ) Fourier transform of a random binary signal 

b T (t) truncated signal of bit) 

b[n] discrete random binary data 

b f [ n] filtered random binary data 

c(x) sound speed 

c 0 reference sound speed 

DFT Discrete Fourier Transform 

d [ n] cosine demodulating signal 

E h average energy per bit 

f c center frequency or carrier frequency 

f s sampling frequency 

f s up-sampling frequency 

FT[] Fourier transform 

FT ' [■] inverse Fourier transform 

g(t - nT h ) rectangular pulse with bit time duration, T b 
G(f ) Fourier transform for a rectangular pulse 


xi 



h(t) ocean impulse response 

h (-t ) time-reversed ocean impulse response 

H(f) ocean frequency response 

H (/) modified ocean frequency response 

h [«] discrete ocean impulse response 

h [«] modified ocean impulse response 

h b [ n] passband ocean impulse response 

H b (/) passband ocean frequency response 

h t \n\ time-reversed ocean impulse response 

h m [n\ mitigated ocean impulse response 

h 0 [ n] baseband ocean impulse response 

H 0 (f ) baseband ocean frequency response 

h p [ n] pre-envelope of passband ocean impulse response 

H (f) pre-envelope of passband ocean frequency response 

h LPF [ n \ FIR lowpass filter 

ISI intersymbol interference 

I interpolation factor 

k [n] cosine modulating signal 

k 0 reference wavenumber 

L number of time and frequency samples used in the transform 

/ index of refraction 


n (t ) additive white Gaussian noise 

n 0 noise component of the output of the correlator 

N bit number of bits in the transmitted binary data 

N 0 /2 two-sided noise power spectral density 

N 0 noise power spectral density 

N noise power 

N sb number of samples per bit 

xii 



P B 

probability of bit error 

P B 

probability of symbol error 

'’b-h) 

conditional probabilities 

pbk) 

conditional probability density function 

Hs,) 

prior probability of the symbols 

pdf 

probability density function 

PSD 

power spectral density 

P b (f) 

PSD for a random binary signal with a rectangular pulse 

PAD 

PSD of BPSK signal 

p[x,t) 

acoustic pressure as a function of position x and time t 

PApskU) 

single-sided power spectral density 

d(0 

received signal 

q[n] 

discrete time received signal 

fiO 

complementary error function 

r 

range 

R 

bit rate (bits/ sec) 

S BPSK 

BPSK signal set 

At) 

BPSK signal 

a(0 

continuous-time real bandpass signal 

s[n\ 

discrete time BPSK signal 

s P ( f ) 

“in-phase” modulation components 

s A T ) 

“quadrature” modulation components 

At) 

complex envelope of BPSK signal 

s(t) 

mitigated signal 

s[n] 

discrete time mitigated signal 

S[k] 

DFT coefficient of s[n] 

S 

signal power 

SNR 

average signal power to average noise power ratio 

xiii 



SNR in decibels 

sound source distribution as a function of position x and time 1 

frequency spectrum of a complex valued low-pass signal 
frequency spectrum of bandpass signal 

At sampling interval 

T h bit time duration or symbol time duration 

t time 

t reduced time 

TL transmission loss 

v ; signal component of the output of the correlator 

w l (7) single waveform for the BPSK signal set 

W bandwidth 

co digital frequency 

x spatial position at which the pressure signal is measured 

y(t ) recovered signal 

y\n\ discrete time recovered signal 

z(T h ) correlator output or decision statistic at bit time duration T b 

z depth 

y 0 decision threshold 

2 

a 0 " white Gaussian noise variance 

<|)(t) sinusoidal carrier 



xiv 



<&(/) 

5(0 

^ (r,z,<p,f) 

x 

9 

0,(0 


Fourier transform of a sinusoidal carrier ty(t) 

delta function 

PE field function 

gross propagation delay 
propagation delay 


azimuthal angle 

phase term of transmitted BPSK signal 


xv 



THIS PAGE INTENTIONALLY LEFT BLANK 


xvi 



ACKNOWLEDGMENTS 


I wish to give recognition to: 

My thesis advisor, Professor Charles W. Therrien, for his time and support in past 
one year. His expertise advice, directions and great help were an enormous support to me 
in achieving this thesis. 

My second reader, Professor Murali Tummala, for his time and help. Especially 
discussing the results with him was a great help to me. 

Professor, Kevin B. Smith, for valuable comments 

My wife, Sug Hyun Park, for her love, devotion and encouragement throughout 
my tour at the Naval Postgraduate School. Thanks to her outstanding skills as a mother 
and a wife. I will always love you, perhaps some days more than others. 

My two sons, Sung Yeop and Sung Moon. Although they are probably too young 
to realize, they constantly kept my daily life in the correct perspective. 

My parents, Tae Young Jung and Cha Soon Cho, for their never-ending love and 
encouragement for me. 

My country and Republic of Korea Navy for sending me to NPS to have this 
education and opportunity. NPS not only gave me a good technical knowledge but also 
awarded me a different way of looking at life. 


xvii 



THIS PAGE INTENTIONALLY LEFT BLANK 



I. INTRODUCTION 


A. GENERAL 

The underwater environment is the most challenging region of the battle space in 
which to communicate effectively. Sound propagation in an underwater communication 
channel becomes a very difficult task due to numerous constraints and limitations 
imposed by the nature of the medium and the environment. The most important 
characteristic s of the underwater environment that put limitations on the performance of 
underwater communication are dispersion caused by the spatial and temperal variability 
and multipath propagation. Multipath propagation produces intersymbol interference 
(ISI) and severe degradation of the bit error performance in digital communication. The 
specific investigation in this thesis is the demodulation of Binary Phase-Shift Keying 
(BPSK) in three different shallow water ocean acoustic channels. 

B. OBJECTIVE 

The objective of this thesis was to study the possibility of using the time-reversed 
ocean impulse response prior to the demodulation of a BPSK signal to mitigate the 
degradation of the bit error performance in multipath propagation. The thesis evaluates 
the bit error performance of the mitigated signal compared to the bit error performance of 
the received signal distorted by the intersymbol interference in three different ocean 
environments in the absence of noise. 

C. ORGANIZATION 

The remainder of this thesis is organized as follows. Chapter II gives the general 
theory related to BPSK binary communication. This describes the definition, waveform 
and constellation of the BPSK signal, the theoretical convergence of the bit error 
probability, and two types of the bit error degradation factors. Chapter III describes the 
experiments and simulation procedure used in this thesis. This includes a time-domain 
and frequency-domain BPSK representation. Power Spectral Density (PSD) for BPSK, a 
description of the Monterey-Miami Parabolic Equation (MMPE) model, the ocean 
impulse response extraction from MMPE model, BPSK demodulation and multipath 
mitigation. Chapter IV provides the main results of the simulation. Included in the results 
are the evaluation for the bit error performance of the BPSK signal, bit error degradation 
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and multipath mitigation for BPSK signal in three different shallow water acoustic 
channels. 
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H. BINARY PHASE-SHIFT KEYING 


This chapter describes the general theory related to binary phase-shift keying 
(BPSK). First, the definition, waveform and its constellation for BPSK modulation are 
explained. Next, as an important criterion of performance in digital communication, the 
probability of bit error is defined and the theoretical convergence of the bit error 
probability is described. Finally, two types of the bit error degradation factors are 
presented. 

A. BINARY PHASE-SHIFT KEYING 

In binary phase-shift keying (BPSK), the phase of a constant-amplitude carrier 
signal is switched between two values according to the two possible “messages” m l and 
m 2 corresponding to binary 1 and 0, respectively. Normally, the two phases are separated 
by 180 degrees. If the sinusoidal carrier has amplitude A c , then the average energy per bit 
1 2 

is B b = — A c T b where T b is the bit time duration, and the transmitted BPSK signal can 
be represented as: 

1 2E, r n 0 < 7 < T, 

J < ( ' >= ,h^ cos [ 2 ’‘/,f+® l (o ■ , „ t 2 - 1 ) 

V T * 1 = \,2 

where the phase term, 0,(0 , will have two discrete values, given by 0 ( =O and 0 2 =n . 
In other words, in BPSK modulation, the modulating data signal shifts the phase of the 
waveform s t ( t) to one of two states, either zero or 71(180°). Instead of using a 

rectangular shaped pulse as shown in Figure 2.1, we can use non-rectangular shaped 
pulses by applying hanning, hamming, and Gaussian windows. However, we shall adhere 
to the simplest case of a rectangular pulse in this thesis. 



Figure 2.1. BPSK Signal in the Time Domain. 
3 








Figure 2.2. BPSK Signal in Frequency Domain. 

Figure 2.1 shows a typical BPSK waveform in the time domain with its abrupt 
phase changes at the symbol transitions. The sketch in Figure 2.2 shows the spectrum of 
the BPSK waveform in frequency domain. For this signal set, the single waveform 

w,(0 = cos(2tt f c t) 0<t<T b (2.2) 


can be defined. Using this signal, the BPSK signal set can be represented as 


S BPSK ~~ { yj^b W l (0 ’ l(0} 


(2.3) 
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B. SIGNAL CONSTELLATION FOR BPSK 


Vector 



M = 2 


S 2 


Si 



Figure 2.3. Signal Constellation for BPSK. 


A BPSK signal waveform, as shown in Figure 2.3, can be geometrically 
represented as vectors or phasors on a polar plot; the vector length corresponds to the 
signal amplitude, and the vector direction corresponds to the phase. Such a representation 
is referred to as the symbol constellation and provides a graphical representation of the 
complex envelope of a BPSK signal for the two possible symbols. The taxis of the 
constellation diagram represents the “in-phase” component of the complex envelope, and 
the Q-axis represents the “quadrature” component of the complex envelope. The distance 
between the signals on a constellation diagram measures separation of the modulation 
waveforms and determines how well a receiver can differentiate between all possible 
symbols when random noise is present. The larger the signal distance, the better the 
chance of correct symbol detection. 

C. PROBABILITY OF BIT ERROR FOR BPSK 

An important measure of performance used for digital modulation is the 
probability of symbol error, P E . It is often convenient to specify system performance by 
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the probability of bit error P B , even when decisions are made on the basis of symbols 
rather than bits. The relationship between P B and P E for orthogonal signaling is: 


P 


B 


P E 

log 2 M 


(2.4) 


where M is the number of symbols. 

For BPSK modulation (M=2), the symbol error probability is equal to the bit error 
probability. When the signals are assumed equally likely and signal y.(t) (i = 1,2) is 
transmitted, the received signal, q(t ), is equal to y (!) + n(t ), where n(t) is modeled as 
additive white Gaussian noise (AWGN). The antipodal signals (signals of equal 
amplitude and opposite polarity) y (!) and s 2 (t) are [cf. Eq. (2.3)]: 


1 (0 = v/EA 17 ! (0 0<t<T b 

(2.5) 

(0 = 0<t<T b 

(2.6) 


The detection and demodulation of the received signal involves a correlation with 
each of the antipodal signals, as explained in Chapter 3. At the end of each symbol 
duration T b , the output of the correlator yields a sample z(T h ), called the detection 
statistic given by 


z(T,) = v,(T,) + h,(TJ i = 1,2 (2,7) 

where v, ( T h ) is the desired signal component and n 0 ( T b ) is the noise component. The 
noise component is a zero mean Gaussian random variable, and thus z (T ;> ) is a Gaussian 
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random variable with a mean of either v, or v 2 depending on whether .v, (f) or s 2 (t) was 
sent. 

The decision stage of the detector will choose the s t (t) with the largest correlator 
output Zj (0, or in this case of equal-energy antipodal signals, the detector, using the 
decision rules, decides: 


Slit) if z(T h ) >y 0 

(2.8) 

s 2 (t) otherwise 

(2.9) 


where y o denotes the decision threshold (equal to 0 for equally probable antipodal 
signals) and z(T h ) is the correlator output at time T b . Two types of detection errors can 
be made. The first type of error takes place if s { (t) is transmitted but the detector 
measures a negative value for z(T fc ) and decides (incorrectly) that signal s 2 (t) was sent. 
The second type of error takes place if signal s 2 (t) is transmitted but the detector 
measures a positive value for z(T 6 ) and decides (again incorrectly) that signal .v, (t) was 
sent. Therefore, the probability of bit error, P B , is the sum: 

F B = P ( S l\ S l) P { S l) + P ( S l \ S l) F { S 2) ( 2 - 10 ) 

where P (y |.v ; j are the conditional probabilities and P (y ) are the prior probability of the 
symbols. For the case when the symbols are equally likely (which is mostly the case): 

p(i,)=p(s 2 ) = | (2.11) 

The expression for the bit error probability then becomes: 
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( 2 . 12 ) 


P B=^ P Uh)+| P ( S I K) 

The conditional probabilities P(s n ly) and P (y |,v 2 ) are found by integrating the 

conditional probability density function (pdf) of the output of the correlator which is 
depicted in Figure 2.4 



Figure 2.4 Conditional Probability Density Functions: p ( z\s ^), p[z |.s\) 

The conditional probability density functions (pdf’s) of the signal with an additive 
white Gaussian noise (AWGN) are: 

i r r _ r 

p(z\s) = —j= —exp - Z ' r V ‘ i = 1,2 (2-13) 
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i.e., the conditional pdfs p (z|.s\) (i = 1,2) are Gaussian with mean value v,, and 

2 2 

variance G 0 where G 0 “ corresponds to the noise variance at the output of the correlator. 
Because of the symmetry of p (z|,s\), the bit error probability of Eq. (2.12) reduces to: 

p B =p ( 5 2h) =p ( 5 ik) ( 2 - 14 ) 

Thus, the probability of a BPSK bit error P B is numerically equal to the area 
under the “tail” of either pdf /((zk) or p(z\s 2 ) that falls on the incorrect side of the 
threshold. We can therefore compute P B by either integrating p (z|.v,) between the limits 
and y o or by integrating /?(z|.v 2 ) between the limits y o and °°, where 
Yo =( v i + v 2 )/ 2 > s the optimum threshold. Hence, 


p b = \p(z\s 2 )dz 

To 


(2.15) 


Using Eq. (2.13), the probability of bit error for BPSK is: 


f 1 

1 avn 

1 

( z - v 2 ^ 

2 

p 

~2 

l a ° J 



dz 


(2.16) 


If we introduce u = 


Z~ V 2 


and o 0 du =dz, the integral simplifies to 


7 1 

r 

du -Q 

( v i ~v 2 "j 

r- — ex P 

(Vi-Vjj/^o '/Ik 


J 


(2.17) 
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where Q{ ) is the complementary error function defined as: 


600 = 



.2 \ 


tiu 


(2.18) 


For equal energy antipodal signaling, such as the BPSK format in Eq. (2.3), the 

receiver output signal components are v l = yjE^ when s l (t) is sent and v 2 = -y[E^ when 

2 

s 2 (t) is sent. For AWGN, the noise variance G 0 of the correlator output can be written 
as N 0 /2, where N 0 /2 is two-sided power spectral density of the noise. Thus, we can 
rewrite probability of bit error as follows: 


Pb = ^ 


-J2E„/N 0 


^/ 2 ^^ 


exp 


.2 3 


Vlu = Q 


2E, 

N n 


(2.19) 


D. BIT ERROR PROBABILITY CONVERGENCE 

In digital communications, we more often use E b /N 0 , a normalized version of 
average signal power to average noise power ratio ( S/N or SNR ), as a figure of merit. 
E b is bit energy and can be described as signal power S times the bit time T b . N 0 is 
noise power spectral density, and can be described as noise power N divided by 
bandwidth W . Since the bit time and bit rate R h are reciprocal, we can replace T b with 
1 jR h and write 


Pa _ ST h _ S/R b (2 20) 

N 0 N/W N/W 

If we use R instead of R b to represent bits/sec, and we rewrite Equation (2.20) 

to emphasize that E b / N 0 is just a version of S/N normalized by bandwidth and bit rate, 
we obtain: 
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( 2 . 21 ) 


E h _ S fW N 
n o ~n{R; 

One of the most important metrics of performance in a digital communication 
system terms is a plot of the bit error probability P B versus E b /N 0 . Figure 2.5 shows the 
“waterfall-like” shape of most such curves. 


P« 



For E b /N 0 > v 0 .P H <P 0 . The dimensionless ratio E h /N 0 is a standard quality 
measure for digital communications system performance. Therefore, the smaller the 
required E b /N 0 , the more efficient is the detection process for a given probability of 
error. 


11 





To optimize (minimize) P B in the context of an AWGN channel and the receiver, 
we need to select the optimum receiving filter and the optimum decision threshold. For 
the binary case, the optimum decision threshold is Y 0 = (v 1 + v 2 )/2 and it was shown in 

Equation (2.17) that this threshold results in P b = 2[(v 1 -v 2 )/2o o ]. Next, for 
minimizing P B , it is necessary to maximize the argument of Q(x) . Thus, we need to 
maximize (v, - v 2 )/2°o or equivalently, maximize 


(v. -v,) , , 

SNR =^ — y 2 - (2.22) 

°o 

where (v, - v 2 ) is the difference of the desired signal components at the filter output at 

time t = T ft , and the square of this difference signal is the instantaneous power of the 

2 2 

difference signal. The SNR can also be maximized by minimizing o 0 ; however, a 0 “ is 
usually not under our control. 

E. TWO TYPES OF BIT ERROR PERFORMANCE DEGRADATION 

The effects of bit error performance degradation can be partitioned into two 
categories. The first is due to a decrease in received signal power or an increase in noise 
or interference power, giving rise to a loss in signal-to-noise ratio or E b /N 0 . The second 
is due to signal distortion, such as might be caused by intersymbol interference (ISI) 
where the tail of a pulse can “smear” into adjacent symbol intervals, thereby interfering 
with the detection process and degrading the error performance. Even in the absence of 
noise, the effects of filtering and channel-induced distortion lead to 1ST Figure 2.6a 
shows the effect of a loss in E b /N 0 ; the solid-line corresponds to the theoretical P B 

versus E h /N 0 curve and the dashed-line corresponds to a degradation effect brought 
about by a loss in E b /N 0 . Figure 2.6b shows an irreducible P B caused by distortion 
(ISI). The solid-line corresponds to the theoretical P B versus E b /N 0 curve and the 
dashed- line corresponds to a degradation effect brought about by ISI instead of a simple 
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loss in signal-to-noise ratio. The bit error rate in this case cannot further be reduced by 
increasing the signal power. 



E h /N 0 (dB) E h /N 0 (dB) 


Figure 2.6. (a) Loss in E b /N 0 . (b) Irreducible P B Caused by Distortion 
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HI. EXPERIMENTAL DESCRIPTION 


In this chapter, we describe the experiments and simulation procedure used in this 
thesis. First, as a part of BPSK signal generation, time-domain and frequency-domain 
BPSK representations and Power Spectral Density (PSD) for BPSK are analytically 
derived. Second, we describe the Monterey-Miami Parabolic Equation (MMPE) acoustic 
simulation model and how to generate the ocean responses in the time-domain and 
frequency-domain from the MMPE model. Third, we describe the ocean environment for 
our simulations and the input data for the MMPE model. Finally, we present the results of 
BPSK demodulation and multi-path mitigation. 

A. BPSK SIGNAL GENERATION 

In Chapter II, we described the properties of BPSK. In this section, time-domain 
and frequency-domain representations for the BPSK waveform are described and related 
to the ocean impulse response, and the power spectral density of the received signal. 
Figure 3.1 shows a diagram for BPSK signal generation. The quantity bit) is a random 

binary signal at baseband. The message data, {a n } represent a set of random variables 


with a n = ±1 and probabilities P(a n - l)= P(a n = -1) = . The function g(t-nT b ) i 


is a 


shaped pulse in general; however we will take it to be a rectangular pulse with bit time 
duration, T h . The signal (j)(t) is a sinusoidal carrier with center frequency f c and 


amplitude A c . The output s(t) is the modulated BPSK signal. 


b(t)= £ a nS{t~nT b ) 

n=—oo * 

5 

i. 

Random binary signal 

BPSK signal 

(^(0 = A,, cos (271./;/) 

Sinusoidal carrier 


Figure 3.1. BPSK Signal Generation. 


15 






The Fourier transform for a rectangular pulse with amplitude one and bit duration 


T b is given by: 


G(f) = 



sintJt/Tj 


(3.1) 


The Power Spectral Density (PSD) for a random binary signal, b(t), will be 
evaluated by first truncating b(t) as follows: 


b T (t)= £ a n g(t-nT h ), -Tj2<t<T 0 /2 

n——N 


(3.2) 


where b T (t) is the truncated signal and T 0 /2 = {N +1/2) T fc . The Fourier transform of the 
truncated signal, b T (t), is then, 


B T (f) = FT[b T (t)]= j^a n FT[g{t-nT b )]= j^a n G{f)e~ mTb 

n=— N n=—N 

= G(/) ■ (3,3) 

n=—N 


where the notation FT[] represents the Fourier transform operation and G (f ) is the 
resulting Fourier transform of g(V). 

The PSD for the random signal bit) can be defined by [Ref 3] 


W)= lim 

T 0 ^°o 
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where B T (f ) represents the Fourier transform of the corresponding truncated signal and 
the overbar represents statistical expectation. Upon substituting (3.3) into (3.4), the PSD 
of b(t ) becomes 


W)= lim 

i n —>00 




a„e 




n——N 


= |G(/)flimf-L|; i )=|G(/)| 2 T lim 


L o n=—N 


2N +1 
T~ 


where we have assumed that the a n are statistically independent. Now, using 
T 0 =2(N + 1/2) T b , we obtain 


W) = |G(/)f ]im 

1 1 N ^> 00 


2N +1 

W+W, 


4-|gc/)| j = t, 


sin(7t/r fc ) 


(35) 


where Eq. (3.l) has been used for G(/). This is the PSD of the signal at ‘baseband.” 

For the modulated BPSK signal we have 

s(t) =b(t)A c cos(2nf c t) (3.6) 


The corresponding PSD is given by 

p s (f) = - A -P,Xf + fc)+ - fc ) (3.7) 

4 4 

(see Ref 3), where P h is the PSD of the baseband signal b(t). Upon substituting (3.5) 
into (3.7) we obtain 
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(3.8) 


P s (f) = 


Ac 2 T b 

4 


\G{f + f c )f+\G{f - f c f 



sin 7t (/ + f c )T 6 

2 

sin tt (/ — f c )T /; 

2 ' 

-< 

4 

Af+fX, J 

+ 

Af-fX \ 



B. MMPE MODEL DESCRIPTION 

The Parabolic Equation (PE) method was introduced into underwater acoustics in 
the early 1970’s by Tappert [Ref 8]. The Monterey-Miami Parabolic Equation (MMPE) 
model [Ref 9] is a numerical model to solve acoustic wave propagation problems in the 
ocean ising the PE method. The MMPE Model is used to compute the response at the 
receiver location of a signal transmitted from a source and traveling through the ocean. 
Figure 3.2 shows the three stages of using this model. In the first stage, we characterize 
the input parameters; in the second, we run the model; and in the last stage, we extract the 
frequency response and ocean impulse response from the outputs of the MMPE model. 


Stage 1: 

^ 

Stage 3: 

Characterize the 

f Stage 2: 

Extract the Ocean 

Input Parameters 

—A Run MMPE ) _^ 

Frequency Response and 

for MMPE 

\. model y 

Ocean Impulse Response: 

model 

' 

H(f) and h(t) 


Figure 3.2. Three Stages to Extract Ocean Response from MMPE Model. 

1. A Brief Description of the MMPE Model 

In order to derive the frequency-domain and time-domain response of the ocean 
underwater communication channel, the general theory behind the parabolic equation 
model will be introduced. Its implementation - the split-step Fourier algorithm for solving 
the PE (PE/SSF) method [Ref 10] is also briefly described. 

In ocean communication, all signals such as p(t) represent acoustic pressure as a 
function of time, and are typically denoted by p[x,t) where x represents the ocean 
coordinates at which the pressure signal is measured. The quantity p[x,t) is often 
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referred to as a space-time signal, since it is a function of the spatial parameters x and 

time 1. When measured at a fixed (receiver) location (say jc = jci) it becomes a function 
of time only. It is transformed into an electrical signal by an acoustic transducer or 
hydrophone. 

The inhomogeneous wave (Helmholtz) equation for the acoustic pressure p(x,t) 
in a medium with sound speed c(x) and density p (xj can be expressed as 



where S (x, / j represents the sound source distribution as a function of both position x 

and time t . In our application (as in most) there is a single source at a particular location 
x=x 0 . Hence s(x,f) is formally a spatial “impulse” at x 0 with time dependence 
s(t) corresponding to the transmitted waveform. 

The parabolic equation model is based on an approximation of the Helmholtz 
wave equation in a cylindrical coordinate system. Because of the ocean’s relative 
shallowness compared to horizontal propagation distance for the majority of 
environments, it is well suited for a description in cylindrical coordinates. 

The signals of interest in underwater communications such as BPSK have 
complicated time dependence; thus analytical solution to the wave equation (3.9) would 
be impossible. However the equation is linear, and a solution can be computed relatively 
easily for sources that have time dependence of the form e j2nfi for a single frequency / . 

Thus for a source with a time dependence of the form Ae j2Kft , the output at a point x = xi 
is dependent on / and can be written as H . (f)e J2nft . This is referred to by people 

xi 

working in ocean acoustics as a “time harmonic” solution to the wave equation The 
quantity H (/) is the ocean frequency reponse and its inverse Fourier transform is the 
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ocean impulse response (Green’s function). The time harmonic solution forms the basis 
of computing the response due to a source with arbitrary time dependence by using 
Fourier analysis/synthesis methods. 

From this point forward the analysis will be carried out in discrete time. This is 
appropriate since numerical techniques (i.e., the DFT) are used to perform the 
computations in both the time and frequency domains. We focus first on the analysis of 
low-pass signals at baseband —W/2< f <W /2 and later extend the discussion to the 
case of bandpass signals such as BPSK centered around some higher carrier frequency. In 
particular, let us sample a signal s(t) at times t= 0, T. 2 T and represent the 

resulting discrete-time signal s[n\= s(nT) in terms of its Discrete Fourier Transform 
(DFT): 


s 


M=tL s M 


.27 ikn 


' k =0 


(3.10) 


Here Spt] is the DFT coefficient representing the frequency component of the signal at 
frequency / = k/LT : 


L-l 


.271 kn 




n=Q 


(3.11) 


and L is the number of time and frequency samples used in the transform. The number of 
samples L and time interval T must of course be chosen to cover the time duration and 
bandwith of the signal. 

With the representation Eq. (3.10) for the source waveform, the output waveform 
at a point x\ can be computed from 
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271 k 


. 27 ikn 

S[k]e~ 


(3.12) 




( 0-7T t- A 


1 i=0 


V ^ J 


where the ocean frequency response has been evaluated at the frequency / = 2nk/L and 
multiplied by the Fourier coefficient S[k ] for the source. 

The ocean frequency response terms that need to be computed by MMPE are 


H 


2iik\ 


k = 0, 1,2, ..., --1 


(3.13) 


which represent frequencies up to the Nyquist “fold-over” frequency f\ /2 = . The 

other terms in the transform represent the “negative” frequencies and may be computed 
from 


f2Kk 3 

= H* 

r 2ii(L-k)' 

l L ) 

X 

L 

\ y 


k 


L 

2 ' 




(3.14) 


Both the forward DFT in Eq. (3.11) needed to compute the frequency coefficients S[k] 
of the source signal and the inverse DFT represented in Eq. (3.12) are performed using 
an FFT program. 

To compute H .(f) we proceed as follows. Assuming a time harmonic solution, 

X v 

the Helmholtz equation takes the form 
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(3.15) 


3p(rip . / ) ^ 
r dr I dr y 


1 d 2 p(r,zW,f) 
- 2 3cp 2 


d 2 p(r z$,f) 

dz 2 


+k 2 l 2 {r,z ,cp)p(r,z,cp,/) = -47UP 0 8(r -r s ), 


where 


PO.Z (p ,/) = pO,zip>? j3IJi 


(3.16) 


The reference wavenumber /c 0 is related to a reference sound speed, c 0 , by 


k 


o 




and the acoustic index of refraction is defined by 


/ = 


c 


0 


c(r,zf\ p) 


(3.17) 


(3.18) 


Note that in this derivation the density variances are neglected. By defining the effective 
index of refraction which contains the appropriate addition terms, the influence of the 
density differences at the water-bottom interface can be included [Ref 8]. 

By assuming that the ocean acts as a waveguide with a cylindrical coordinate 
system, acoustic energy is mainly propagated outward from a source in the horizontal 
direction. Therefore, the pressure field can be approximated by 


p(r,z,cp,/) = 'T(r,z,cp,/)H 0 (1) {k Q r) 


(3.19) 
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where 'P (r, z.cp,/) is a slowly varying envelope function, or PE field function and 
H o <n (V) is the zero-th order Hankel function of the first kind of the outgoing acoustic 
wave. Taking advantage of the far-field (k 0 r » l) asymptotic approximation of the 
Hankel function, Eq. (3.19) can be rewritten as 


p(r ,z,cp,/) 


^P(r,z,cp,/>^ 


S 


(3.20) 


normalized such that at the reference range r = R 0 , \p\ = P 0 . Substituting Eq. (3.20) into 
Eq. (3.15) and dropping the source term on the right hand side gives 


3 2v P 1 3 2v P 3 2x P 

- —h ]2k () -1 -— H--—h 

dr 2 dr r 2 3cp" dz 2 


(i 2 -i)+-L 

V ’ 4r 2 


¥=0 


(3.21) 


Neglecting the azimuthal coupling and the far-field terms, and dropping the first term due 
to the slow modulation of the envelope function, Eq. (3.2l) can be rewritten as 


j_dV 
k 0 dr 


1 3 2v P 

2 k 2 dz 2 


j (‘ 2 -!)*■ 


(3.22) 


Now, by defining the operators 


1 

T =- 

op 2k, 


2 dz 2 


V 


(3.23) 


and 
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(3.24) 


U 


op 



Eq. (3.22) can be rewritten as 


j_<W_ 

k 0 dr 


= (T +U W . 

V op °P I 


(3.25) 


Equation (3.25) is known as the “standard” parabolic equation (SPE) [Ref 8], with 

accurate solutions limited to a half width of ±15° for the propagation angle and 
represents the complete description of the outgoing acoustic energy propagation in the 
waveguide. In order to extend this limit to ± 40°, a higher order wide-angle parabolic 
equation (WAPE) approximation [Ref 11] is used. Its operators are defined by 


T'wAPE — 1 Jl + 


kjd? 


(3.26) 


and 


U\VAPE l) ■ 


(3.27) 


The MMPE uses the split-step Fourier (SSF) method in order to solve the 
parabolic equation numerically. This algorithm integrates the solution in range by 
applying the T WAPE and the U WAPE operators in the wavenumber (k 7 )-domain and the 
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spatial (z)-domain, respectively, where each operator is a scalar multiplier. In the k_- 


A 

space, the wide angle T wape operator to simplify calculation is defined as 


A 

T WAPE 



(3.28) 


The field function at range r + Ar is expressed as 


T^r + Ar, 


• ft 


Ar. 
~jk n - 1 


]ka—UwAPE( r ' z ) 


'F(az) 


(3.29) 


Since the output data of the MMPE model is the field function in the form of 
equation (3.29), we substitute Eq. (3.29) into Eq. (3.20). The ocean frequency response 
(with respect to pressure) at range r and depth z is then given by 


H [rz) (f)=p(r, Z ,f) = -j='¥{r, Z ,f) e ik ° r 


(3.30) 


As described earlier, for broadband signals, the output results require knowing 
Hi .)(/) for all frequencies in the bandwidth of the signal [see Eq. (3.12)]. Thus the 

MMPE model must be run for all frequencies within the band. The model produces 
frequency domain results which are written to an output file. A post processing program 
performs the synthesis specified in Eq. (3.12). 

2. Signal Representation in MMPE 

The time intervals encountered in ocean acoustics can be large. In particular the 
time due to propagation delay from source to receiver driving which no sound is heard at 
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the receiver can be very large compared to the duration of the signal when it finally 
arrives (see Figure 3.3) 



Received 

Signal 



Figure 3.3. Pulse Arriving after Propagation Delay X . 

To obtain the most accurate results, the analysis time needs to be focused in a small time 
window surrounding the received signal as shown in Figure 3.3. 
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The signals of interest are also often restricted in frequency as well as time. For 
example, a BPSK transmission has a frequency spectrum similar to that shown in Figure 
3.4. 



Such signals are referred to as bandpass signals; the spectral power is assumed negligible 
outside of an interval ±W /2 around the center frequency or carrier frequency f c . Again, 
for accuracy and reduced computation in processing, it is most appropriate to focus in a 
frequency window (frequency band) where the spectrum is non-zero. 

MMPE deals with the signal processing in an efficient way by permitting the user 
to specify a time window and frequency band for the processing. The formulation for 
such processing is described below. 

a. Time Windowing 

The time windowing in MMPE is relatively simple to explain. The 

7 271 /— 

frequency response H (r ^(f) of Eq. (3.30) contains a term e jk ° r =e * which 
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corresponds to the gross propagation delay x = r/c 0 from source to receiver in the ocean. 
MMPE drops this term, so that the modified frequency response is given by 




(3.31) 


The time domain signals resulting from this modified frequency response correspond to a 
reduced time scale t defined as 


■ r 

t -t-x-t - 


(3.32) 


and referred to as reduced time . MMPE uses this reduced time for its output results. 
Absolute time (relative to the source) can be obtained by adding back in the gross 
propagation delay x =r/c 0 . 

Without the modification to the frequency response specified in Eq. 
(3.31), some time response would be produced within the analysis window of the DFT. 
However, this response would appear as a circularly shifted version of the true correct 
response and the time scale would not be correct. This effect is known as “time aliasing.” 
b. Processing of Bandpass Signals 

A continuous-time real bandpass signal such as that depicted in Figure 3.4 
can be written in the general form 


s b (0 = S P {t)cos2K f c t-s q (t) sin2jt f c t 
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Here s (t) and s (t) are known as the “in-phase” and “quadrature” modulation 
components. The bandpass signal can be equivalently expressed as 


s b (t)= Re s(t)e j2Kfc ' 


(3.33) 


where 


5 (0 = 5 p {t) + js g (t) 


is known as the complex envelope . The complex envelope is a complex-valued low-pass 
signal with spectrum as depicted in Figure 3.5. 



Figure 3.5 Spectrum of a Complex Valued Lowpass Signal Corresponding to 

Bandpass Signal of Figure 3.4. 
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The complex envelope may be obtained from the original bandpass signal by discarding 
the negative-frequency part of the spectrum, moving the positive-frequency part of the 
spectrum from a center point of / = f c to / = 0 (“baseband”) and scaling by a factor of 
2 . 

The advantage of using complex envelope representation for the bandpass 
signal is that the processing can be focused in the bandpass region by processing the 
complex envelope at baseband. Let us denote the sampled version of the complex 

envelope as 


s[n] = s(nT) = s p ( nT) + js q ( nT) 


(The notation s\n\ was previously used to denote a sampled real band-limited baseband 
signal; the only difference now is that s[n] is complex.) The complex baseband signal 
can be represented as in Eq. (3.10) where S[k] are the frequency components at 
baseband defined by Eq. (3.1l). However a frequency in the DFT of 271 k/L represents 
an actual frequency in the ocean of 2n k/L + f c . Hence the ocean frequency response 
H- (/) must be computed by MMPE at this higher frequency (27 lk/L+ f c ). The time- 
sampled complex envelope of the received signal is then given by 


q[n\ 


1 L -1 

-yH 

Lh x 


ik , l r ., 

— + fc S l k ] e 


.2.% kn 


(3.34) 


(Compare this to Eq. (3.12)). The sampled real bandpass signal at the receiver is thus 
given by 


30 



q h [/■?] = 2Re|p[//]e ;2,t/ ‘" r J 


(3.35) 


3. Input Parameters for the MMPE Model 

The input data for the MMPE model fall into three main categories: the main 
input file, the source data, and the environmental data. Each of these is described below. 

The main input file is called pefiles . inp: it is the only filename that must 
exist. All other input files are specified in this file. Figure 3.6 shows that the first several 
entries of this main file simply define the other input files. Line 1 specifies the source 
data file. Lines 2 to 6 specify files that provide the environmental data, while line 7 
specifies the output data file (a binary file). Following these data files, the user specifies 
the number of points in depth, minimum depth, maximum depth (at line 8), the number of 
points in range, minimum range, maximum range (at line 9). The final line in this file 
represents the vertical FFT size, the range step, the maximum depth of the calculation, 
and the reference sound speed. 


Environmental 

data 


0 i Q § % i g * n M /► | I i I i ® I i i | Stack: f - X 



1 

Bane of source data input file :pesrc.inp 


p 

Name of ssp data input file :pessp.inp 


3 

Name of bottom bathy data file rpebath.inp 

-< 

4 

Name of bottom properties file :pebotprop.inp 

w 

5 

Name of deep bottom bathy file ipedbath.inp 


u 

Name of deep bottom props file :pedbotprop.inp 


7 

Name of output data file :posfc400bw200c.bin 


8 

nzout, depmin[m], depmax[m] :128 0 200 


9 

nrout, rnpin[km], rnpax[km] : 1000 0 10.0 


10 

nz, dr[km], depcalc[m], c0[m/s] :256 0.0 0.0 1500.0 


11 




il J±l 


Figure 3.6. pefiles . inp File of the Main Input File. 


The source data input file depicted in Figure 3.7 defines all of the source 

information. There are two source types which are available: a wide-angle source which 

approximates a point source, and a vertical line array which approximates a continuous 

line array. For this work, we will use the wide-angle point source approximation The 
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first line in Figure 3.7 represents the source cbpth. The second line denotes the array 
length; an array length of zero indicates a point source. For a wide-angle source, the D/E 
angle is ignored since a point source cannot be steered. The remaining lines represent the 
center frequency, the bandwidth, and the number of the frequencies at which the ocean 
response is to be evaluated (size of the DFT). This number must be a power of two. 



Figure 3.7. pesrc . inp File of the Source Data. 

The environmental data consists of five input files. The pessp . i np file contains 
the sound speed profile data and is illustrated in Figure 3.8. As shown there, the first line 
contains two numbers indicating the number of azimuthal radials and the total azimuthal 
aperture. The second fine contains a single number indicating the number of sound speed 
profiles. The following line denotes the range of the current profile and the number of 
sound speed values. Finally, the profile is defined by pairs of the depth and sound speed. 

The other environmental files are similar in structure. Their detailed descriptions 
are contained in the Ocean Acoustic Library web page 
http://oalib.saic.com/PE/mmpe.html which is supported by the U.S. Office of Naval 
Research. The pebath. inp file contains the bathymetry of the water/bottom interface; 
the pebotprop. inp file contains the acoustics parameters of the medium just below 
the water/bottom interface; the pedbath . inp file defines the “deep” layer bathymetry 
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beneath the water/sediment interface; and finally, the pedbotprop. inp file contains 
the acoustic properties of the deep layer. 



The last file, the output data file is a binary file produced by MMPE. It is read by the 

postprocessing program to produce final results 

C. OCEAN ENVIRONMENTAL CHARACTERIZATION 

For our simulation studies of underwater communication, we focused on shallow 
water environments. Three different ocean environments which have three different 
sound speed profiles (SSP) were simulated to test the performance of underwater 
communication using a BPSK signal. The complete set of MMPE input files for all three 
cases are given in Appendix A. 

1. Case 1: Positive SSP Gradient 

For this case, the source is located at a range of 0 m and at a depth of 5 m. The 
sound speed is range-independent with a positive, linear gradient SSP of 1497 m/s at the 
surface and 1499 m/s at a depth of 100m. The density of water is assumed to be 1.0 

a 

g/cm . The bottom is chosen to be a ‘bumpy’ bottom with the following depths: 102m at 
0 km; 101 m at 2.5 km; 99m at 5 km; 102 m at 7.5 km; and 99 m at 10 km. Its 

compressional sound speed is 1700 m/s, the sound speed gradient 1 s~ l , density 1.5 
g/cm 3 , and compressional attenuation 0.1 dB/km/Hz. Figure 3.9 shows the positive 
sound speed profile versus depth and the corresponding sound transmission loss plot 
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(versus range and depth) for a source at depth 5 m and frequency of 400 Hz. 
Transmission loss is defined as 


TL = 201og[p(l, z s )/p(r,z )] 


(3.36) 


where p(r,z ) is the acoustic pressure amplitude and p(\,zj is the reference acoustic 
pressure amplitude measured at range 1 m (from the source) and at the source depth z s . 

The color in Figure 3.9 represents the transmission loss in dB according to the scale 
shown on the right. 




Figure 3.9. (a) Positive SSP (b) Sound Propagation at Frequency 400 Hz, Source 5 m. 


2. Case 2: Strong Negative SSP Gradient 

For this case, the source is located at a range of 0 m and at a depth of 30 m. The 
sound speed is range-independent with a strong negative, bilinear gradient (downward 
refraction) sound speed profile of 1528 m/s at the surface and 1510 m/s at a depth of 
50m, and 1489 m/s at the depth of 100 m. The density of water is assumed to be 1.0 
g/cm. The bottom is chosen to be a ‘bumpy’ bottom with the same depths as in the 
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previous case. Its compressional sound speed is 1700 m/s, the sound speed gradient 1 
.v~ 1 , density 1.5 g/cm 3 , and compressional attenuation 0.1 dB/km/Hz. Figure 3.10 shows 
the strong negative sound speed profile versus depth and the sound transmission loss for a 
source at depth 30 m, and frequency of 400 Hz. 



GtMid 5pe»d. |rrV5| 



Figure 3.10.(a) Negative SSP (b) Sound Propagation at Frequency 400 Hz, Source 30 m. 


3. Case 3: Negative SSP Gradient Below Surface Duct 

For this case, the source is again located at a range of 0 m and at a depth of 5 m. 
The sound speed is range-independent with a negative, bi-linear gradient sound speed 
profile with surface duct of 1492 m/s at the surface and 1500 m/s at a depth of 40m, and 
1489 m/s at the depth of 100 m. The density of water is assumed to be 1.0 g/cm. The 
bottom is chosen to be a ‘bumpy’ bottom with the same depths as in the previous cases. 

Its compressional sound speed is 1700 m/s, the sound speed gradient 1 s' 1 , density 1.5 
g/cm, and compressional attenuation 0.1 dB/km/Hz. Figure 3.11 shows the positive 
sound speed within the surface duct, the negative sound speed profile below the surface 
duct, and the sound transmission loss for a source at depth 5 m and frequency of 400 Hz. 
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Figure 3.11. (a) SSP (b) Sound Propagation at Frequency 400 Hz, Source 5 m. 

D. BPSK DEMODULATION AND DETECTION 

At the receiver location, the received signal waveforms arrive distorted due to 
noise and intersymbol interference in the underwater communication channel. We need to 
demodulate the received signal in order to recover the transmitted signal and recover the 
binary data. Figure 3.12 shows two basic steps in the demodulation and detection process. 
Step 1, the waveform-to-sample transformation, consists of the demodulator followed by 
a sampler. At the end of each symbol time duration T h , the sampler produces an output 

z( T fc ) which in the absence of noise is proportional to the energy of the received symbol. 
Step 2 is a decision process where z(T 6 ) is compared to a threshold y 0 to decide if the 
received data represents a binary 1 or binary 0. 

1. Correlation Receiver 

As mentioned in Chapter II, there are two main types of degradation factors in the 
performance of underwater communication, namely noise and dispersion introduced by 
the ocean channel with its multipath environment. In this section, we ignore the 
degrading factor produced by the ocean impulse response and assume that the only 
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performance degradation is due to AWGN with zero mean. In next section, the effect of 
the ocean impulse response will be added to the received signal. 



Figure 3.12. Two Basic Steps in the Demodulation/Detection of the Received Signal. 

In the process of demodulation, the received signal is reduced to a single random 
variable z(T fo ) sometimes called a “detection statistic”. In the absence of distortion due 
to the ocean, the received signal can be expressed as: 


q{t) = s(t -X) + n(t) 0 <t <T 


(3.37) 


where q(t) represents the received signal, s(t-x) is the transmitted delayed BPSK 
signal, and n(t) is AWGN with zero mean and variance, o 0 2 . 

As shown in Figure 3.13, the recovered signal is formed by multiplying the input 
signal by two local sinusoidal carriers .v, (t) and s 2 (t) , assumed synchronized with the 
received signal to obtain 


y l (t) = q(t)-s i (t)= q(t} cos(2x f c t) (3.38) 
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and 


y 2 (0 = q(t)-s 2 (t) = -q(t)-cos(2nf c t) (3.39) 



Figure 3.13. Correlator Receiver. 


The recovered signals (t ) and y 2 ( t ) are integrated over the bit interval to 


obtain 


z,-( T *) = j yS t ) dt = \ q(t)$(t)dt 1=1,2 (3.40) 


and the outputs of the integrations are subtracted to form the detection statistic 


z(T*)= (T& ) — ^2 ( T *) 


(3.41) 
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In the decision stage, z(T b ) is compared to an optimum threshold as follows 


z(T b ) * ^=Y 0 (3-42) 

where Vj is the signal component of z(T ; J when .v, (t) is transmitted, and v 2 is the signal 
component of z(T h ) when s 2 (t) is transmitted. The threshold level y 0 defined by 
To = ( v i + v 2 )/ 2 is the optimum threshold for minimizing the probability of bit error [Ref 
1]. For equal energy, equally likely antipodal signals, where s, (t) = —s 2 (t) and v , = -v 2 , 
the decision rule becomes 


z(J b ) J 0 =0 (3.43a) 

Assuming that s^t) corresponds to binary 1 and s 2 (t) is binary 0, the decision 
rule thus reduces to 

decide binary 1 if Zj(T fc ) > z 2 (T 6 ) (3.436) 

decide binary 0 otherwise 

E. MULTIPATH MITIGATION 

In this section, we assume the received signal is further distorted by intersymbol 
interference due at least in part to multipath propagation in the channel. We need to use 
some multipath mitigation technique to compensate for the degraded signal. One such 
technique is time reversal. In this method the received signal is convolved with the time- 
reversed impulse response of the ocean before applying it to the demodulation. As shown 
in Figure 3.14, the received signal has the following form 

q(t) - s(t ) * hit) + n(t) 
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(3.44) 



where h{t) is the ocean impulse response in the time-domain. 


Local sinusoidal carrier 


s l (t)-s 2 (t) = 2cos( 2nf c t) 



Received signal Recovered signal 


A A 

q(t) = s(t) * h(t) + n(t) y(t) = s(t)-2cos(2n f c t), where s(t) = s(t)* h(t)*h(-t) 


Figure 3.14 Multipath Mitigation for the Distorted Signal in Underwater Environment. 

The distortion in the signal (represented by the convolution s(t)*h(t)) is 

mitigated by applying the time-reversed ocean impulse response to the received signal. In 
the absence of noise, the resulting mitigated signal is 

s(t) = q{t)*h (~t)= s(t)*h(t)*h (—t) (3.45) 

The convolution of h{t) with h ( -t ) produces an impulse-like signal. Hence the 
result of Eq. (3.45) is a tendency to restore the distorted signal to its original condition. 

The recovered signal is then formed by multiplying the mitigated signal by the two local 
sinusoidal carriers s { (t) and s 2 (t) to obtain 
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(3.46) 


y(t)=s(t)-(s,(t)-s 2 (t)) = s(t)-2cos(2nf,t) 


The recovered signals y(t) are integrated over the bit time duration to form the 
detection statistic. This is compared to a threshold as before to produce the recovered 
binary data. 
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IV. SIMULATION RESULTS 


A. EVALUATION OF BIT ERROR PROBABILITY FOR BPSK SIGNAL 

In this section, the parameters (bandwidth, sampling frequency, bit rate, samples 
per bit, interpolation factor and up-sampling frequency) for the BPSK signal and a Finite 
Impulse Response (FIR) filter are defined, and the effect of the AWGN on the bit error 
performance is evaluated. 

1. Evaluation for BPSK Parameters 

An important theorem of communication is based on the assumption of a strictly 
bandlimited channel, i.e. one in which no signal power whatever is allowed outside the 
band of interest. For our work, we need to define the bandwidth for the BPSK signal 
transmission. The single-sided power spectral density for the BPSK signal (also known as 
the analytic signal) is given by 


1 BPSK 


(f) 


At T, 


sintt {f~f c )T h 

*(/-/,) T* 


(4.1) 


This follows from Eq. (3.8) by dropping the terms for negative frequencies. This 

power spectral density is depicted in Figure 4.1 and is seen to consist of a main lobe and 
smaller sidelobes. Although there are many criteria for measuring bandwidth, for our 
digital communication, we are constrained to two bandwidth criteria namely Null-to-null 
bandwidth and Power bandwidth. 

The Null -to-null bandwidth is given by W = 2R = 2/T b where R is the 

position of the first null relative to the center frequency (see Figure 4.1). The sampling 
frequency corresponding to this definition of bandwidth is given by 

f s = 2R = 2/T, (4.2) 
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Thus, the bit rate for this bandwidth definition is R = W/2 = \/T h and the number 
of samples per bit is N sh = f : T b = 2 . 



Figure 4.1. The Power Spectral Density for BPSK Signal. 


The Power bandwidth defines the frequency band in which 99% of the total 
power resides. This bandwidth has been adopted by the Federal Communication 
Commission (FCC Rules and Regulations Section 2.202) and states that the occupied 
bandwidth is the band that leaves exactly 0.5% of the signal power above the upper band 
limit and exactly 0.5% of the signal power below the lower band limit. Thus for this 
definition 99% of the signal power is inside the band. For the BPSK signal this 
bandwidth is given by W = 20.56 R = 20.56/T,, . The sampling frequency for the power 
bandwidth is given by 


f s = 20.56R = 20.56/T, 
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(4.3) 



























Thus, the bit rate for this bandwidth is R = W/ 20.56 =1 /T b and the number of 
samples per bit is N sb = f s T h = 20.56 . 

Since the carrier frequency is generally much higher than the bandwidth of the 
baseband signal, a random binary signal must be interpolated to reconstruct a bandlimited 
waveform without error. The interpolation factor is determined as 


I = ceil 


2 (/ r + W) 

f, 


(4.4) 


where I represents the interpolation factor and ceil rounds towards plus infinity. The 
increased sampling rate is then 


f',=V, (4.5) 

where f s represents the up-sampling frequency. 

The interpolated random binary signal must then be filtered to remove the 
unwanted spectral energy above the band. A digital Finite Impulse Response (FIR) filter 
is used. The FIR filter is designed using a Hamming window with the following 
parameters. 


1.0(27t IT/2) 

|G) „ =-;-= passband digital frequency 

f f 

2.0(271 IT/2) . 1J? 

|co. =-- — t -- = stopband digital frequency 


/, 


N,=- 


871 


0 ) - 0 ) 

s p 


■ filter order (number of coefficients is N h +1) 
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The filter parameters are shown for a general FIR filter designed using the 
Hamming window method in Figure 4.2. The variable co , is the digital frequency defined 
as 


oo = 


2 nf 

fs 


(4.6) 



Figure 4.2. General FIR Filter (Hamming Window Design). 
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In the following experiments, the frequency response of the BPSK signal for the 
Null-to-null bandwidth and Power bandwidth is shown. Figure 4.3 shows the frequency 
response for 1000 bits of a BPSK signal which is sampled according to the Null-to-null 
bandwidth and has parameters R =100 bits/sec, f c =400 Hz, / 1 =200 Hz. N sb = 2, 1=6 

and f s =1200 Hz. Figure 4.4 shows the frequency response for 1000 bits of a BPSK 
signal which is sampled according to the Power bandwidth criterion with parameters 
R=100 bits/sec, f =4000 Hz, /,=: 2100 Hz, 7=6, f' s = 12600 Hz and N sb =2\ 
= ceil (20.56). The complete procedures to generate the BPSK signal are explained in 
Appendix B. 



Figure 4.3. Frequency Response for BPSK Signal Sampled by Null-to-null 

Bandwidth. 



As shown in Figure 4.3 and Figure 4.4, sampling according to the Power 
bandwidth is more reliable than the Null-to-null bandwidth since 99% of the signal power 
is within the bandwidth. However, since the most energy is contained within the Null-to- 
null bandwidth, it is sufficient to work with the Null-to-null bandwidth for our work. In 
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the next subsection, the bit error performance under AWGN for BPSK signal having 
Null-to-null bandwidth is evaluated. 

2. Influence of AWGN on Bit Error Probability 


For evaluation of the bit error performance under AWGN, the bit error probability 
was tested for various values of the average signal power to average noise power ratio 
( SNR ). The average signal power was kept constant. 

E 

The average signal power to average noise power ratio is defined as SNR =— 

a A/ 


where E,=^ e ^ 


b - is the average signal power, o 0 is the white Gaussian noise variance, 


and At = ——. The SNR is thus given by 

N.„ 


SNR = 


a, 2 t. 


2ct 0 A t 


(4.7) 


or in decibels as 


' a 2 T ^ 

SNR (d B)=101og ( SNR ) = 10 log c 


(4.8) 


The theoretical value of the bit error probability P B from Eq. (2.19) is: 


P b=Q 


f [2eT A 


JK J 


q{4snr) = q 


f ' A 2 T,, ^ 


2a At 

v ’ 0 j 


(4.9) 
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In the following experiment, the parameters for the simulation to evaluate the bit 
error performance under AWGN for a BPSK signal are as follows: N bit =10000 bits, 

R =100 bits/sec, W =200 Hz, f c =400 Hz, f s =200 Hz, N sb =2. The noise variance a 0 “ is 
varied from 0.01 to 2 in increments of 0.1 and from 3 to 60 in increments of 2. 


BER vs SNR for BPSK : Nbit=10000, Nsb=2, fs=200 Hz, fc= 400 Hz 



Figure 4.5. BER versus SNR (dB) . 


The experimental bit error rate estimates for various values of SNR for the Null- 
to-null bandwidth is shown in Figure 4.5. It can be seen that the experimental values 
follow the theoretical values given by Eq. (4.9) but are over all slightly lower. This is 

due to the effect of the FIR filter which removes some of the noise. Since the effect of 
AWGN on BER for BPSK signal is consistant with the theoretical values, our concern in 
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ocean communication will be now with the degradation due to ocean impulse response 
and the multipath than the degradation due to the additive noise. 

B. BIT ERROR DEGADATION AND MULTIPATH MITIGATION FOR 
BPSK SIGNAL IN A SHALLOW WATER ENVIRONMENT 

In this final section, the bit error performance degradation due to the ocean 
impulse response is evaluated. This impulse response is generated from the MMPE 
model in the three different ocean environments as described in Chapter III. In all three of 
these cases, we used the same bandwidth (200 Hz) and center frequency (400 Hz) for the 
BPSK signal. In all cases the performance is compared for two situations. First, the 
received signal is applied directly to the demodulation with no prior mitigation steps and 
the BER is evaluated as a function of range and depth from the transmitter. Secondly, the 
received signal is convolved with the time-reversed ocean impulse response before 
applying it to the demodulation. This mitigation step tends to compress the spreading 
caused by the ocean impulse response. The BER is then evaluated and compared to that 
of the unmitigated situation. 

1. Bit Error Performance Results for a Positive SSP Gradient 

The parameters of this environment are given in Chapter III (see Figure 3.9). 
First, by choosing an ocean impulse response extracted from the MMPE model at a 
chosen depth of 50 m and range 5 km, the influence of the ocean impulse response on the 
bit error rate at this chosen depth and range is evaluated. Then the time-reversed ocean 
impulse response is convolved with the received signal to compensate for the distorted 
signal and the bit error rate is evaluated again. In the following, we describe this one 
specific case in detail. Later, in section B the evaluation is performed at many chosen 
depths and ranges of this ocean environment. 

a. Investigation Of The One Specific Case 

The parameters of the generated random binary data are as follows: 
N bit =20 bits, R = 100 bits/sec, N sb = 24, f s =2400 Hz This discrete random binary data 

b[n\ is filtered to remove an unwanted signal (i.e., above bandwidth, W =200 Hz) by a 

FIR filter. The filtered random binary data is b j [n] = b [«] *h LPF [«] (see Appendix B). 

This filtered random binary data b f [ n] is modulated by cosine modulating signal, 
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cos 


r 2n f n ^ 

-— , having a carrier frequency of f c =400 Hz and the sampling interval 

V J 

Al = 1//' . This entire procedure simulates the sampling of an analog BPSK waveform 
(which is what would actually be present in the water) at 2400 samples/sec. The 
modulated BPSK signal, as shown in Figure 4.6 in time and frequency domains, is 
2 nf c n 


>[n] = b f [n]• cos 


fs 



The modified ocean frequency response H (/) extracted from the 

MMPE model at a chosen depth of 50 m and range 5 km corresponds to a center 
frequency of f c = 400 Hz, a bandwidth of W =300 Hz and a sampling frequency of 

f s = 300 Hz. To obtain the ocean impulse response at the same sampling rate and 
bandwidth as the BPSK signal, the ocean impulse response from MMPE is up-sampled 


by a factor of 8, 


I = 


2(3x/,) 


= 8 , and filtered to remove an unwanted signal (i.e., above 


bandwidth, 200 Hz) by a FIR filter (see Appendix C). Figure 4.7 shows the passband 
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ocean frequency response H b (f ) (magnitude and phase) and ocean impulse response 
h b [n ]. The complete procedure to generate the passband ocean impulse response at this 
high sampling frequency is explained in Appendix C. 
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Figure 4.7. Passband Ocean Frequency Response (Magnitude and Phase) and Ocean 

Impulse Response. 


The BPSK signal s[n] is convolved with the ocean impulse response 
h b [n] by linear convolution in the time domain. The received signal is 
q[n] = v[n] */j 6 [h]. This simulates the distortion caused by propagation through the 
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ocean. Figure 4.8 shows the distorted received signal in time and frequency domains at 
the chosen receiver location. 


The recovered signal is obtained by multiplying the received signal q[n] 


by the cosine modulating signal, cos 
2 nf c n 


r 2xf c n A 
f 

V J s / 


Thus, the recovered signal is 


y[/z] =<?[«]• cos 


fs 


. This recovered signal is filtered to remove unwanted 


frequency components above the bandwidth by a FIR filter and integrated over the bit 
time duration to form the detection statistic. This is compared to a threshold to produce 
the recovered binary data. Figure 4.9 shows that the recovered binary data has a very high 
bit error rate (55 %) due to the ocean impulse response when compared to the transmitted 
binary data. 



Figure 4.8. Received Signal in Time Domain and Frequency Domain. 
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Figure 4.9. Comparison of Recovered Binary Data with Transmitted Binary Data for 

the Received Signal (55 % error rate). 


To mitigate this high bit error rate due to the ocean impulse response, the 
time- reversed ocean impulse response is applied to the received signal. The time- 
reversed ocean impulse h t [n\ is obtained by reversing the passband ocean impulse 

response h b [n\ in time (i.e. h t [«] = h b [—«]). Figure 4.10 shows the passband time- 

reversed ocean frequency response (magnitude and phase) and ocean impulse response. 
The phase is the negative of the phase in Figure 4.7. 

To show the mitigated ocean impulse response, the time-reversed ocean 
impulse response is convolved with the ocean impulse response. The mitigated ocean 
impulse response is thus given by h m = h b [ n] *h r [«] = h b [ n\ * h b [ —n]. Figure 4.11 shows 
that in the absence of noise the mitigated ocean impulse produces an impulse-like signal. 
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Figure 4.10. Time Reversed Ocean Frequency Response (Magnitude and Phase) and 

Ocean Impulse Response. 
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Figure 4.11. Mitigated Ocean Impulse Response and Frequency Response. 


The mitigated signal is obtained by convolving h b [ -n ] with the received 
signal to obtain s [ n\ = q [ n] * h b [-»]. The recovered signal is obtained by multiplying the 

mitigated signal ,v[»] by cosine to obtain y [n] = s [n] • cos 

same procedures to recover the transmitted binary data described before. The recovered 
binary data is compared to the transmitted binary data to compute BER. Figure 4.12 
shows the mitigated signal in time and frequency domains obtained by convolving the 
time-reversed ocean impulse response with the received signal. Figure 4.13 shows that 
the mitigated procedure was able to reduce the bit error rate from 0.55 % to 0 %. 


2 nf c n 

fs 


and following the 
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Figure 4.12. Mitigated Signal in Time and Frequency Domain. 



Figure 4.13. Comparison of Recovered Binary Data with Transmitted Binary Data for 

Mitigated Signal Showing Zero Error. 
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b. Results For Many Chosen Ranges and Depths 

In the following experiment, to obtain more complete results of simulation 


in this ocean environment, the simulation was performed for the received signal and the 
mitigated signal to evaluate the bit error performance for many chosen depths and ranges. 
The parameters for this simulation are the same as in the previous section except the 
number of bits was increased from 20 to 10,000 (N hlt =10,000). The ocean frequency 

responses H (/) were extracted from the MMPE model at depths of 5 m to 95 m in 
increments of 5 m and ranges of 0.5 km to 9.5 km in increments of 0.5 km with a 
sampling frequency f s = 300 Hz, bandwidth W =300 Hz and center frequency f c =400 
Hz. The experimental procedures followed were as described before. 

As shown in Figure 4.14, the recovered binary data set for the received 
signals at chosen depths and ranges as described above has a high bit error rate 
(average=0.4946) due to the distortion of the transmitted BPSK signal due to the 
multipath propagation resulting from the positive SSP gradient. Figure 4.15 however 
shows that the recovered binary data set for the mitigated signals obtained by convolving 
the time-reversed ocean impulse response with the received signal has a low bit error rate 
(average=0.0451) compared to the bit error rate of the unmitigated received signals. The 
results of this simulation at many chosen depths and ranges is consistent with the result of 
simulation at the former chosen depth of 50 m and range 5 km. We conclude therefore 
that by using the time-reversed ocean impulse response, we can compensate for the 
degradation in the bit error performance due to the multipath propagation in this 
simulated ocean environment. 
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Figure 4.14. Bit Error Performance for Received Signal at Chosen Depths 5 to 95 m 
and Ranges 0.5 to 9.5 km for a Positive SSP Gradient. 
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Figure 4.15. Bit Error Performance for Mitigated Signal at Chosen Depths 5 to 95 m 
and Ranges 0.5 to 9.5 km for a Positive SSP Gradient. 


2. Bit Error Performance Results for a Strong Negative SSP Gradient 

In the following experiment, for a shallow water acoustic channel of having a 
strong negative SSP gradient, the simulation was performed for the received signal and 
the mitigated signal to evaluate the bit error performance for many chosen depths and 
ranges. The SSP and TL plots are shown in Figure 3.10a and 3.10b. The signal 
parameters and grid of range and depth values were the same as for the previous case and 
the experimental procedures followed were as described before. As shown in Figure 4.16, 
the recovered binary data set for the received signals at chosen depths and ranges as 
described above has a high bit error rate (average=0.4981) due to the distortion of the ISI 
resulting from the strong negative sound speed profile. Figure 4.17 shows that the 
recovered binary data set for the mitigated signals however has a low bit error rate 
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(average=0.0415) compared to the bit error rate of the unmitigated received signals. We 
conclude therefore that by using the time-reversed ocean impulse response, we can 
compensate for the degradation in the bit error performance due to the multipath 
propagation in this ocean environment as well. 



Figure 4.16. Bit Error Performance for Received Signal at Chosen Depths 5 to 95 m 
and Ranges 0.5 to 9.5 km for a Strong Negative SSP Gradient. 
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Figure 4.17. Bit Error Performance for Mitigated Signal at Chosen Depths 5 to 95 m 
and Ranges 0.5 to 9.5 km for a Strong Negative SSP Gradient. 

3. Bit Error Performance Results for a Negative SSP Gradient below 
Surface Duct 

In the following experiment, for a shallow water acoustic channel having a 
negative SSP gradient below surface duct, the simulation was performed for the received 
signal and the mitigated signal to evaluate the bit error performance for many chosen 
depths and ranges. The SSP and TL plots for this environment are shown in Figure 3.11a 
and 3.11b. The signal parameters and grid of range and depth values were the same as in 
the previous case and the experimental procedures followed were as described before. As 
shown in Figure 4.18, the recovered binary without mitigation has a high bit error rate 
(average=0.4989) due to the multipath and the resulting 1ST Figure 4.19 shows that the 
recovered binary data set for the mitigated signals however has a low bit error rate 
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(average=0.0398) compared to the bit error rate of the unmitigated received signals. We 
conclude that for this environment as well, by using the time-reversed ocean impulse 
response, we can compensate for the degradation in the bit error performance due to the 
multipath propagation. 
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Figure 4.18. Bit Error Performance for Received Signal at Chosen Depths 5 to 95 m 
and Ranges 0.5 to 9.5 km for a Negative SSP Gradient below Surface Duct. 
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Figure 4.19. Bit Error Performance for Mitigated Signal at Chosen Depths 5 to 95 m 
and Ranges 0.5 to 9.5 km for a Negative SSP Gradient below Surface Duct. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The primary goal of this thesis was to investigate mitigating the degradation on 
the bit error performance of the BPSK signal by convolving the time-reversed ocean 
impulse response with the received signal distorted by the multipath propagation in three 
different shallow water acoustic channels. The simulation was validated by a comparison 
of the results of the bit error performance from the received signal and those from the 
mitigated signal. In our experiments, it was possible to reduce the high bit error rate as 
follows. For a positive SSP gradient, the average bit error rate decreased from 0.4946 to 
0.0451. For a strong negative SSP gradient, the average bit error rate decreased from 
0.4981 to 0.0415. For a negative SSP gradient below the surface duct, the average bit 
error rate decreased from 0.4989 to 0.0398. Thus, it was possible to improve the bit error 
performance of a BPSK signal by using the time-reversed ocean impulse response. 

The experiments assumed that the correct ocean impulse response was used at 
each position of the receiver. The sensitivity to range and depth or incorrect ocean 
impulse response was not investigated in any quantitative manner. There were some 
indications, however, of a fair amount of sensitivity to change in position. 

B. RECOMMENDATIONS 

Although our experiment was not able to produce a sufficiently low bit error rate 
needed to achieve an effective underwater communication in the ocean environment, the 
algorithm may be of some use if appropriate error correction coding is employed to 
reduce the bit error rate. Experiments with a time-reversed ocean impulse response 
showed reasonable success in reducing the bit error rate of the distorted BPSK signal due 
to the intersymbol interference. The combination of appropriate other filtering and 
coding/decoding with this time-reversed ocean impulse response may therefore be worthy 
of further investigation. 


65 



THIS PAGE INTENTIONALLY LEFT BLANK 


66 



APPENDIX A. MMPE INPUT FILES FOR THREE DIFFERENT 
OCEAN ENVIRONMENTAL CASES 


This appendix gives the complete set of MMPE input files for three different 
ocean environmental cases as described in Chapter III. The MMPE input files of the 
environmental data except the input file of the sound speed profile and pesrc . inp file 
of the source data are same for three cases. 

A. MMPE INPUT FILES FOR POSITIVE SSP GRADIENT 
1. pef iles . inp File of the Main Input File 


D 

fe « r> | m /► | @ @ I ^ m a m | ««* 


i 

Name of source data input file 

:pesrc.inp 

▲ 

2 

Name of ssp data input file 

:pessp.inp 


3 

Name of bottom bathy data file 

:pebath.inp 


4 

Name of bottom properties file 

:pebotprop.inp 


5 

Name of deep bottom bathy file 

:pedbath.inp 


6 

Name of deep bottom props file 

:pedbotprop.inp 


7 

Name of output data file 

:posfc400bw300b.bin 


8 

nzout, depmin[m], depmax[m] 

:10 0 200 


9 

nrout, rngmin[km], rngmaxfkm] 

o 

o 
■—1 

o 

o 

H 


10 

nz, dr[km], depcalc[m], cO[m/s] 

:256 0.0 0.0 1500.0 


11 





<1 


M 


2. pesrc. inp File of the Source Data 
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3 


pessp. inp File of the Environmental Data 


D 


& & 

i 

1 0 . 


2 

i 


3 

0. 101 


4 

0 

1497 

5 

1 

1497 

6 

2 

1497 

7 

3 

1497.1 

8 

4 

1497.1 

9 

S 

1497.1 

10 

6 

1497.1 

11 

7 

1497.1 

12 

8 

1497.2 

13 

9 

1497.2 

14 

10 

1497.2 

15 

11 

1497.2 

16 

12 

1497.2 

17 

13 

1497.3 

18 

14 

1497.3 

19 

IS 

1497.3 

20 

16 

1497.3 

21 

17 

1497.3 

22 

18 

1497.4 

23 

19 

1497.4 

24 

20 

1497.4 

25 

21 

1497.4 




o c* 


/► ® £) , €) 


Stack: . 2£j 




4. pebath. inp File of the Environmental Data 
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5. 


pebotprop. inp File of the Environmental Data 



6. pedbath. inp File of the Environmental Data 



7. pef iles . inp File of the Environmental Data 
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B. MMPE INPUT FILES FOR STRONG NEGATIVE SSP GRADIENT 
1. pef iles . inp File of the Main Input File 


D h m 

<& % d 10 c* 

m /► 

% 

^ Ur 41 m 

Stack: | 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Name of source data input file 
Name of ssp data input file 
Name of bottom bathy data file 
Name of bottom properties file 
Name of deep bottom bathy file 
Name of deep bottom props file 
Name of output data file 
nz out, depmin[m], depmax[m] 
nr out, rngmin[km], rngmax[km] 
nz, dr[km], depcalc[m], c0[m/s] 

<1 


:pesrc.inp 
:pessp.inp 
:pebath.inp 
:pebotprop.inp 
:pedbath.inp 
:pedbotprop.inp 
:negfc400bw300a.bin 
:128 0 200 
:1000 0 10.0 
:256 0.0 0.0 1500.0 



2. pessp. inp File of the Environmental Data 


□ 

& s * 

S X 

1 

1 0. 


2 

1 


3 

0. 101 


4 


o is; 

5 

1 

1527.6 

6 

2 

1527.3 

7 

3 

1526.9 

8 

4 

1526.6 

9 

5 

1526.2 

10 

6 

1525.8 

11 

7 

1525.5 

12 

8 

1525.1 

13 

9 

1524.8 

14 

10 

1524.4 

15 

11 

1524 

16 

12 

1523.7 

17 

13 

1523.3 

18 

14 

1523 

19 

15 

1522.6 

20 

16 

1522.2 

21 

17 

1521.9 

22 

18 

1521.5 

23 

19 

1521.2 

24 

20 

1520.8 

25 

21 

1520.4 


<1 


d ^ c* 


/► €) €) I €) 


Stack: 


.r 1 
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C. MMPE INPUT FILES FOR NEGATIVE SSP GRADIENT BELOW 
SURFACE DUCT 

1. pef iles . inp File of the Main Input File 


□ g* B # 1 1 rn H, « o. | A /> I ® £| | ® ®r id id 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Name of source data input file 
Name of ssp data input file 
Name of bottom bathy data file 
Name of bottom properties file 
Name of deep bottom bathy file 
Name of deep bottom props file 
Name of output data file 
nzout, depmin[m], depmax[m] 
nrout, rngmin[km], rngmax[km] 
nz, dr[km], depcalc[m], cO[m/s] 

<1 


:pesrc.inp 
:pessp.inp 
:pebath.inp 
:pebotprop.inp 
:pedbath.inp 
rpedbotprop.inp 
:surfc400bw300a.bin 
:128 0 200 
:1000 0 10.0 
:256 0.0 0.0 1500.0 


Stack: [sa^T X| 


i >r 


2. pessp. inp File of the Environmental Data 



1 SB S 

(#4 /► 

@ © 

^ if iH *£) 

Stack: Basi X_ 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


1 0. 


▲ 

1 

0. 101 

0 1528 


1 

1527.6 


2 

1527.3 


3 

1526.9 

— 

4 

1526.6 


5 

1526.2 


6 

1525.8 


7 

1525.5 


8 

1525.1 


9 

1524.8 


10 

1524.4 


11 

1524 


12 

1523.7 


13 

1523.3 


14 

1523 


15 

1522.6 


16 

1522.2 


17 

1521.9 


18 

1521.5 


19 

1521.2 


20 

1520.8 


21 

1520.4 

▼ 
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APPENDIX B. DETAILED SIGNAL PROCESSING STEPS 


This appendix gives the complete procedures from the generation of the BPSK 
signal to the demodulation of the received signal for an ideal digital communication 
channel (i.e., without any distortion) 

A. RANDOM BINARY DATA GENERATION 

The generated random binary data b[n\ has parameters N bit =20 bits, R =100 

bit/sec, N sb = 24, f s = 2400 Hz and T h =0.01 sec. Figure B-l shows the random binary 
data in time and frequency domains. 



Freq Response of Binary Data 



Figure B-1. Random Binary Data in Time Domain and Frequency Domain. 
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B. FILTERING 


This random binary data is filtered to remove an unwanted spectral energy above 
bandwidth {W = 200 Hz). The FIR filter h LPF [n\ is designed using a Hamming window 

as described in chapter 4. The filtered binary data is b f [n] =b[n\ *h LPF \n\ . Figure B-2 
shows the frequency response of FIR filter h LPF \n\ and the filtered binary data b f \n]. 



Figure B-2. Frequency Response of the FIR Filter and the Filtered Binary Data. 


C. GENERATION OF BPSK SIGNAL 


The filtered binary data b f [n] is modulated by a cosine modulating signal, 


k [n] = cos 


2 nf c n 

fs 


having a carrier frequency of f c = 400 Hz. The modulated BPSK 


signal is ,s [n] = b j [n ] ■ k [n]. Figure B-3 shows the modulated BPSK signal in the time 
and frequency domains. 
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Figure B-3. BPSK Signal ,y[/z] in Time Domain and Frequency Domain. 


D. DEMODULATION OF BPSK SIGNAL, FILTERING 

The transmitted BPSK signal s[n] is demodulated by multiplying by a cosine 

2 nf c n 


demodulating signal, d \n\ = cos 


fs 


The demodulated signal is given by 


y[n\ =s[n\ d [n] and produces a double frequency component centered at ±800 Hz as 
shown in Figure B-4. The signal y [ n ] is filtered again to remove unwanted spectral 
energy outside bandwidth by using same FIR filter described above. The demodulated 
and filtered signal is y f [n] = y[n]*h LPF [n \. Figure B-4 shows the frequency response of 

the demodulated BPSK signal y[/z] before lowpass filtering. Figure B-5 shows the 
demodulated and filtered signal y f [n] in time and frequency domains. 
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Figure B-4. Demodulated BPSK Signal y [ n ] in Frequency Domain. 



Frsquancy Posponso oi DomoduiJtod and Lv^assad BPSK Signal 



Figure B-5. Demodulated and Filtered BPSK Signal in Time and Frequency Domain. 


E. BER COUNTING 

This demodulated and filtered signal y f \n\ is integrated over the bit time 
duration to form the detection statistic and then compared to an optimum threshold as 
described in Eq. (3.43a) of Chapter III and by this decision rule, we decide binary 1 or 0. 
Thus, the recovered binary data is compared to the transmitted binary data. Figure B-6 
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shows that the bit error rate of the recovered binary data for an ideal digital 
communication signal goes to zero. 



Figure B-6. Comparison between Recovered Binary Data and Transmitted Binary 

Data. 
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APPENDIX C. COMPLETE PROCEDURES TO GENERATE THE 
PASSBAND OCEAN IMPULSE RESPONSE AND FREQUENCY 
RESPONSE FROM THE MM PE MODEL 


This appendix gives the complete procedures to generate the passband ocean 
impulse response and frequency response at a high sampling frequency from the modified 

ocean frequency response H (/) extracted from MMPE model. 

A. MODIFIED OCEAN FREQUENCY RESPONSE AND OCEAN IMPULSE 
RESPONSE FROM MMPE MODEL 

The modified ocean frequency response H (/) extracted from MMPE model as 

mentioned in Eq. (3.31) of Chapter III has IT = 300 Hz, f c = 400 Hz, f s =300 Hz 

and n . =256, where n f represents the number of frequency components or FFT size. 

Let us consider a chosen depth of 50 m and range 5 km in the ocean environment having 
a positive SSP gradient. Figure C-l shows the modified ocean frequency response 

H (/) and ocean impulse response h [//] obtained by taking the inverse DFT of 


B. BASEBAND OCEAN IMPULSE RESPONSE BY PADDING ZEROS 

To obtain the ocean impulse response at the same sampling frequency as the 
BPSK signal, we need to pad zeros between lower part and upper part of the modified 

ocean frequency response H (/). In this case, to obtain 8 times sampling frequency 

(f s , =8 f s ) of the original sampling frequency of f s =300 Hz, we need to pad 


VP IT 

zeros between -3 f c and 3 f c , except for the interval — —< /< —. These zeros should 


be equally spaced as A/ 


IT 


n f-\ 


■ 1.1765 Hz. Figure C-2 shows the baseband ocean 


frequency response H 0 (f ) and ocean impulse response h o [n]. 
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Figure C-l. 


Modified Ocean Frequency Response H (/) and Impulse Response 




Figure C-2. 


Baseband Ocean Frequency Response H a (f ) and Impulse Response 
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C. FILTERING AND PRE-ENVELOPE OF THE OCEAN IMPULSE 
RESPONSE 

Since the baseband ocean response W = 300 Hz, we need to filter it to have the 
same bandwidth as the BPSK signal (200 Hz). The FIR filter h LPF [n] is designed using a 

Hamming window as described in chapter 4. Figure C-3 shows the frequency response of 
the FIR filter and the filtered baseband ocean frequency response. 

The filtered baseband ocean impulse response h f o [ n] is multiplied by 
exp(-j27t f c t) to obtain the pre-envelope of the passband ocean impulse response, 
h p [ n] = h j o [n ] ■ exp ( —j2n ft ). Figure C-4 shows the frequency response of the pre¬ 
envelope of the ocean impulse response H (f). 


Fi»quency Response of LFf with Pessbao!l= 10D Hz. Siopl>en<f= 330 Hz 



■ 150-1-1-1-1-L 


o xo 400 boo aoo lino 1230 

Frequency (Hz) 



Figure C-3. Frequency Response of FIR Filter and Filtered Baseband Ocean 

Frequency Response. 
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Figure C-4. Pre-Envelope Ocean Frequency Response II p (/). 


D. PASSBAND OCEAN IMPULSE RESPONSE 

Since the pre-envelope of the passband ocean impulse response is 

h [n] = h b [n] + jh b [n\ where h b [n\ is the passband ocean impulse response and /?, [n] is 
the Hilbert transform of h b [n] (see Ref 13). The passband ocean impulse response is 
obtained by taking the real part of h p [n] and by scaling by a factor of 2. The passband 
ocean impulse response is defined as h b [n] = 2Re [h p [n]). The passband frequency 

response H b (/) is obtained by doing the Fourier transform of h b \n\. Figure C-5 shows 
the passband ocean frequency response. 
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